{smcl}
{com}{sf}{ul off}{txt}{.-}
      name:  {res}<unnamed>
       {txt}log:  {res}/N/project/suicide_study/pnas_replication/results/log/mi_4.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res}20 Aug 2020, 05:05:12
{txt}
{com}. 
. if ("`model'" == "logit"){c -(}
.         use "${c -(}home_dir{c )-}/data/processed/suicide_reg_v1_raw.dta", clear
. {c )-}
{txt}
{com}. 
. if ("`model'" == "mi") {c -(}
.         use "${c -(}home_dir{c )-}/data/processed/suicide_reg_v1_imputed_M10.dta", clear  
. {c )-}
{txt}
{com}. 
. * for now, we use the following simple survey weights 
. if ("`model'" == "mi") {c -(}
.         mi svyset `geo_type' [pw=ObsWgt0] 
{res}
      {txt}pweight:{col 16}{res}ObsWgt0
          {txt}VCE:{col 16}{res}linearized
  {txt}Single unit:{col 16}{res}missing
     {txt}Strata 1:{col 16}<one>
         SU 1:{col 16}{res}county
        {txt}FPC 1:{col 16}<zero>
{p2colreset}{...}
{res}{com}. {c )-}
{txt}
{com}. else {c -(}
.         svyset `geo_type' [pw=ObsWgt0]  
. {c )-}
{txt}
{com}. 
. * margins for each category
. program margin_interact 
{txt}  1{com}.         args X Y k model
{txt}  2{com}.         sum `X', d 
{txt}  3{com}.         local gap = (`r(max)' - `r(min)') / `k' 
{txt}  4{com}.         if ("`model'" == "logit") {c -(}
{txt}  5{com}.                 margin `Y', at(`X' = (`r(min)' (`gap') `r(max)')) predict(pr)
{txt}  6{com}.         {c )-} 
{txt}  7{com}.         else if ("`model'" == "mi") {c -(}
{txt}  8{com}.                 mimrgns `Y', at(`X' = (`r(min)' (`gap') `r(max)')) predict(pr)
{txt}  9{com}.         {c )-}       
{txt} 10{com}. end 
{txt}
{com}. 
. program mchange_mi
{txt}  1{com}.         args X k model
{txt}  2{com}.         if ("`model'" == "logit") {c -(}
{txt}  3{com}. 
.                 if ("`k'" == "continuous") {c -(}
{txt}  4{com}.                         sum `X' if e(sample), d 
{txt}  5{com}.                         margin, at(`X' = (`r(min)' `r(max)')) post predict(pr)
{txt}  6{com}.                         mlincom  2 - 1, decimal(7) stat(all)    
{txt}  7{com}.                 {c )-} 
{txt}  8{com}.                 else if ("`k'" == "binary") {c -(}
{txt}  9{com}.                         sum `X' if e(sample), d 
{txt} 10{com}.                         margin, at(`X' = (`r(min)' `r(max)')) post predict(pr)
{txt} 11{com}.                         mlincom  2 - 1, decimal(7) stat(all)
{txt} 12{com}.                 {c )-} 
{txt} 13{com}.                 else if ("`k'" == "categorical") {c -(}
{txt} 14{com}.                         margin `X' if e(sample)==1 , at() pwcompare predict(pr) 
{txt} 15{com}.                 {c )-}
{txt} 16{com}.         {c )-}
{txt} 17{com}. 
.         else if ("`model'" == "mi") {c -(}
{txt} 18{com}. 
.                 if ("`k'" == "continuous") {c -(}
{txt} 19{com}.                         sum `X' , d 
{txt} 20{com}.                         mimrgns, at(`X' = (`r(min)' `r(max)')) post predict(pr)
{txt} 21{com}.                         mlincom  2 - 1, decimal(7) stat(all)
{txt} 22{com}.                 {c )-} 
{txt} 23{com}.                 else if ("`k'" == "binary") {c -(}
{txt} 24{com}.                         sum `X' , d 
{txt} 25{com}.                         mimrgns, at(`X' = (`r(min)' `r(max)')) post predict(pr)
{txt} 26{com}.                         mlincom  2 - 1, decimal(7) stat(all)
{txt} 27{com}.                 {c )-} 
{txt} 28{com}.                 else if ("`k'" == "categorical") {c -(}
{txt} 29{com}.                         mimrgns `X' , at()   pwcompare predict(pr)      
{txt} 30{com}.                 {c )-}
{txt} 31{com}.         {c )-}
{txt} 32{com}. end
{txt}
{com}. 
. * create some dummy codings
. tab Race5, gen(race5_nh)

      {txt}Race5 {c |}      Freq.     Percent        Cum.
{hline 12}{c +}{hline 35}
          1 {c |}{res}107,831,845       79.53       79.53
{txt}          2 {c |}{res} 14,631,144       10.79       90.32
{txt}          3 {c |}{res}  2,867,304        2.11       92.44
{txt}          4 {c |}{res}  2,753,311        2.03       94.47
{txt}          5 {c |}{res}  7,499,129        5.53      100.00
{txt}{hline 12}{c +}{hline 35}
      Total {c |}{res}135,582,733      100.00
{txt}
{com}.         rename race5_nh1 White_nh 
{res}{txt}
{com}.         rename race5_nh2 Black_nh 
{res}{txt}
{com}.         rename race5_nh3 AIAN_nh 
{res}{txt}
{com}.         rename race5_nh4 AsPI_nh 
{res}{txt}
{com}.         rename race5_nh5 Hispanic
{res}{txt}
{com}. 
. tab MarStat5, gen(ms)

   {txt}MarStat5 {c |}      Freq.     Percent        Cum.
{hline 12}{c +}{hline 35}
          1 {c |}{res} 75,689,936       55.83       55.83
{txt}          2 {c |}{res} 10,166,664        7.50       63.32
{txt}          3 {c |}{res} 13,925,509       10.27       73.60
{txt}          4 {c |}{res}  2,609,308        1.92       75.52
{txt}          5 {c |}{res} 33,190,855       24.48      100.00
{txt}{hline 12}{c +}{hline 35}
      Total {c |}{res}135,582,272      100.00
{txt}
{com}.         rename ms1 Marrd5
{res}{txt}
{com}.         rename ms2 Widow5
{res}{txt}
{com}.         rename ms3 Divor5
{res}{txt}
{com}.         rename ms4 Separ5
{res}{txt}
{com}.         rename ms5 NvMar5
{res}{txt}
{com}. 
. tab AgeGrp4, gen(ag) 

    {txt}AgeGrp4 {c |}      Freq.     Percent        Cum.
{hline 12}{c +}{hline 35}
          1 {c |}{res} 21,231,562       15.66       15.66
{txt}          2 {c |}{res} 39,043,312       28.80       44.46
{txt}          3 {c |}{res} 47,351,260       34.92       79.38
{txt}          4 {c |}{res} 27,956,599       20.62      100.00
{txt}{hline 12}{c +}{hline 35}
      Total {c |}{res}135,582,733      100.00
{txt}
{com}.         rename ag1 Age_15_24
{res}{txt}
{com}.         rename ag2 Age_25_44
{res}{txt}
{com}.         rename ag3 Age_45_64
{res}{txt}
{com}.         rename ag4 Age_65_Up
{res}{txt}
{com}. 
. destring St, replace 
{txt}St: all characters numeric; {res}replaced {txt}as {res}byte
{txt}
{com}. 
. * set-up equations
. local religion Rat_GC_ProE Rat_GC_ProM Rat_GC_ProB Rat_GC_Cath Rat_GC_Jew Rat_GC_Oth
{txt}
{com}. local contextual_control Rat_Poverty Rat_Mig_Cum Pop_Den
{txt}
{com}. 
. local religion Rat_GC_ProE Rat_GC_ProM Rat_GC_ProB Rat_GC_Jew Rat_GC_Oth
{txt}
{com}. local contextual_control Rat_Poverty Rat_Mig_Cum Pop_Den
{txt}
{com}. 
. local demographics_raw i.Female c.RAT_Female i.AgeGrp4 c.RAT_AgeGrp4_2 c.RAT_AgeGrp4_3 c.RAT_AgeGrp4_4 i.Race5 c.RAT_Race5_2 c.RAT_Race5_3 c.RAT_Race5_4 c.RAT_Race5_5 i.BornUSA c.RAT_BornUSA i.MarStat5 c.RAT_MarStat5_2 c.RAT_MarStat5_3 c.RAT_MarStat5_4 c.RAT_MarStat5_5
{txt}
{com}. local demographics_same i.Female c.std_same_prop_Sex i.AgeGrp4 c.std_same_prop_AgeGrp4 i.Race5 c.std_same_prop_Race5 i.BornUSA c.std_same_prop_BornUSA i.MarStat5 c.std_same_prop_MarStat5 
{txt}
{com}. local demographics_inter i.Female##c.std_same_prop_Sex i.AgeGrp4##c.std_same_prop_AgeGrp4 i.Race5##c.std_same_prop_Race5 i.BornUSA##c.std_same_prop_BornUSA i.MarStat5##c.std_same_prop_MarStat5
{txt}
{com}. 
. 
. if (`model_version' == 1){c -(}
.         local model_eq i.Year `demographics_raw' `contextual_control' `religion' 
.         local margin_demographics Female RAT_Female AgeGrp4 RAT_AgeGrp4_2 RAT_AgeGrp4_3 RAT_AgeGrp4_4 Race5 RAT_Race5_2 RAT_Race5_3 RAT_Race5_4 RAT_Race5_5 BornUSA RAT_BornUSA MarStat5 RAT_MarStat5_2 RAT_MarStat5_3 RAT_MarStat5_4 RAT_MarStat5_5 
. {c )-}
{txt}
{com}. if (`model_version' == 2){c -(}
.         local model_eq i.Year `demographics_raw' `contextual_control' `religion' i.UnEmpl c.RAT_UnEmpl i.PhysProb c.RAT_PhysProb
.         local margin_demographics Female RAT_Female AgeGrp4 RAT_AgeGrp4_2 RAT_AgeGrp4_3 RAT_AgeGrp4_4 Race5 RAT_Race5_2 RAT_Race5_3 RAT_Race5_4 RAT_Race5_5 BornUSA RAT_BornUSA MarStat5 RAT_MarStat5_2 RAT_MarStat5_3 RAT_MarStat5_4 RAT_MarStat5_5 UnEmpl RAT_UnEmpl PhysProb RAT_PhysProb
. {c )-}
{txt}
{com}. if (`model_version' == 3){c -(}
.         local model_eq i.Year `demographics_same' `contextual_control' `religion' 
.         local margin_demographics Female std_same_prop_Sex AgeGrp4 std_same_prop_AgeGrp4 Race5 std_same_prop_Race5 BornUSA std_same_prop_BornUSA MarStat5 std_same_prop_MarStat5 
. {c )-}
{txt}
{com}. if (`model_version' == 4){c -(}
.         local model_eq i.Year `demographics_same' `contextual_control' `religion' i.UnEmpl c.std_same_prop_UnEmpl i.PhysProb c.std_same_prop_PhysProb
.         local margin_demographics Female std_same_prop_Sex AgeGrp4 std_same_prop_AgeGrp4 Race5 std_same_prop_Race5 BornUSA std_same_prop_BornUSA MarStat5 std_same_prop_MarStat5 UnEmpl std_same_prop_UnEmpl PhysProb std_same_prop_PhysProb
. {c )-}
{txt}
{com}. if (`model_version' == 5){c -(}
.         local model_eq i.Year `demographics_inter' `contextual_control' `religion' 
. {c )-}
{txt}
{com}. if (`model_version' == 6){c -(}
.         local model_eq i.Year `demographics_inter' `contextual_control' `religion' i.UnEmpl##c.std_same_prop_UnEmpl i.PhysProb##c.std_same_prop_PhysProb
. {c )-}       
{txt}
{com}. * test how long it would take.
. * mi estimate: svy: mean Suic 
. * local demographics i.Female c.RAT_Female i.AgeGrp4 c.RAT_AgeGrp4_2 c.RAT_AgeGrp4_3 c.RAT_AgeGrp4_4 i.Race5 c.RAT_Race5_2 c.RAT_Race5_3 c.RAT_Race5_4 c.RAT_Race5_5 i.BornUSA c.RAT_BornUSA i.MarStat5 c.RAT_MarStat5_2 c.RAT_MarStat5_3 c.RAT_MarStat5_4 c.RAT_MarStat5_5
. * mi estimate: svy: logit Suic i.St `demographics' UnEmpl RAT_UnEmpl PhysProb RAT_PhysProb
. 
. * main effects : margins
. if ("`model'" == "mi"){c -(}
. 
.         mi estimate: svy: logit Suic i.St `model_eq', or 
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Survey: Logistic regression{col 47}Number of obs{col 65}= {res}  11,814,307

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918
{txt}{col 47}Average RVI{col 65}= {res}      0.0188
{txt}{col 47}Largest FMI{col 65}= {res}      0.3307
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}       79.03
{txt}{col 47}        avg{col 65}= {res}      865.90
{txt}{col 47}        max{col 65}= {res}      914.77
{txt}Model F test:{ralign 16: {res:Equal FMI}}{col 47}F({res}  51{txt},{res}  914.3{txt}){col 65}= {res}      838.31
{txt}Within VCE type: {ralign 12:{res:Linearized}}{col 47}Prob > F{col 65}= {res}      0.0000

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}        Suic{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 10}St {c |}
{space 10}8  {c |}{col 14}{res}{space 2}-.0042586{col 26}{space 2}  .092176{col 37}{space 1}   -0.05{col 46}{space 3}0.963{col 54}{space 4}-.1851599{col 67}{space 3} .1766426
{txt}{space 9}13  {c |}{col 14}{res}{space 2}-.3876931{col 26}{space 2} .0928146{col 37}{space 1}   -4.18{col 46}{space 3}0.000{col 54}{space 4}-.5698475{col 67}{space 3}-.2055387
{txt}{space 9}21  {c |}{col 14}{res}{space 2}-.5399436{col 26}{space 2} .0971564{col 37}{space 1}   -5.56{col 46}{space 3}0.000{col 54}{space 4}-.7306191{col 67}{space 3} -.349268
{txt}{space 9}24  {c |}{col 14}{res}{space 2}-.5659662{col 26}{space 2} .0936554{col 37}{space 1}   -6.04{col 46}{space 3}0.000{col 54}{space 4}-.7497712{col 67}{space 3}-.3821613
{txt}{space 9}25  {c |}{col 14}{res}{space 2}-.8857383{col 26}{space 2}  .093791{col 37}{space 1}   -9.44{col 46}{space 3}0.000{col 54}{space 4}-1.069809{col 67}{space 3}-.7016675
{txt}{space 9}34  {c |}{col 14}{res}{space 2} -.796781{col 26}{space 2} .0947736{col 37}{space 1}   -8.41{col 46}{space 3}0.000{col 54}{space 4}-.9827801{col 67}{space 3}-.6107819
{txt}{space 9}35  {c |}{col 14}{res}{space 2} .1560838{col 26}{space 2} .0972092{col 37}{space 1}    1.61{col 46}{space 3}0.109{col 54}{space 4}-.0346957{col 67}{space 3} .3468632
{txt}{space 9}37  {c |}{col 14}{res}{space 2}-.3591295{col 26}{space 2} .0947141{col 37}{space 1}   -3.79{col 46}{space 3}0.000{col 54}{space 4}-.5450116{col 67}{space 3}-.1732473
{txt}{space 9}40  {c |}{col 14}{res}{space 2}-.5008854{col 26}{space 2} .1029172{col 37}{space 1}   -4.87{col 46}{space 3}0.000{col 54}{space 4}-.7028677{col 67}{space 3} -.298903
{txt}{space 9}41  {c |}{col 14}{res}{space 2}-.2807343{col 26}{space 2} .0936794{col 37}{space 1}   -3.00{col 46}{space 3}0.003{col 54}{space 4}-.4645859{col 67}{space 3}-.0968828
{txt}{space 9}44  {c |}{col 14}{res}{space 2} -.730422{col 26}{space 2} .0946151{col 37}{space 1}   -7.72{col 46}{space 3}0.000{col 54}{space 4}-.9161101{col 67}{space 3}-.5447338
{txt}{space 9}45  {c |}{col 14}{res}{space 2}-.3878547{col 26}{space 2} .0957385{col 37}{space 1}   -4.05{col 46}{space 3}0.000{col 54}{space 4}-.5757474{col 67}{space 3}-.1999619
{txt}{space 9}49  {c |}{col 14}{res}{space 2} .4912192{col 26}{space 2} .1585365{col 37}{space 1}    3.10{col 46}{space 3}0.002{col 54}{space 4} .1800818{col 67}{space 3} .8023567
{txt}{space 9}51  {c |}{col 14}{res}{space 2}-.2324933{col 26}{space 2} .0919496{col 37}{space 1}   -2.53{col 46}{space 3}0.012{col 54}{space 4}-.4129514{col 67}{space 3}-.0520353
{txt}{space 9}55  {c |}{col 14}{res}{space 2}-.4997238{col 26}{space 2} .0966312{col 37}{space 1}   -5.17{col 46}{space 3}0.000{col 54}{space 4}-.6893688{col 67}{space 3}-.3100788
{txt}{space 12} {c |}
{space 8}Year {c |}
{space 7}2006  {c |}{col 14}{res}{space 2}-.0364229{col 26}{space 2} .0183438{col 37}{space 1}   -1.99{col 46}{space 3}0.047{col 54}{space 4}-.0724238{col 67}{space 3}-.0004219
{txt}{space 7}2007  {c |}{col 14}{res}{space 2} .0377063{col 26}{space 2} .0187117{col 37}{space 1}    2.02{col 46}{space 3}0.044{col 54}{space 4} .0009834{col 67}{space 3} .0744292
{txt}{space 7}2008  {c |}{col 14}{res}{space 2} .0475266{col 26}{space 2} .0203233{col 37}{space 1}    2.34{col 46}{space 3}0.020{col 54}{space 4} .0076407{col 67}{space 3} .0874125
{txt}{space 7}2009  {c |}{col 14}{res}{space 2}-.0062114{col 26}{space 2} .0196069{col 37}{space 1}   -0.32{col 46}{space 3}0.751{col 54}{space 4}-.0446921{col 67}{space 3} .0322694
{txt}{space 7}2010  {c |}{col 14}{res}{space 2}-.0412063{col 26}{space 2}  .021474{col 37}{space 1}   -1.92{col 46}{space 3}0.055{col 54}{space 4}-.0833521{col 67}{space 3} .0009395
{txt}{space 7}2011  {c |}{col 14}{res}{space 2} .0163382{col 26}{space 2} .0230416{col 37}{space 1}    0.71{col 46}{space 3}0.478{col 54}{space 4}-.0288831{col 67}{space 3} .0615595
{txt}{space 12} {c |}
{space 4}1.Female {c |}{col 14}{res}{space 2}-1.335525{col 26}{space 2} .0155227{col 37}{space 1}  -86.04{col 46}{space 3}0.000{col 54}{space 4} -1.36599{col 67}{space 3}-1.305061
{txt}std_same_p~x {c |}{col 14}{res}{space 2}-.0156892{col 26}{space 2}  .007495{col 37}{space 1}   -2.09{col 46}{space 3}0.037{col 54}{space 4}-.0303985{col 67}{space 3}-.0009798
{txt}{space 12} {c |}
{space 5}AgeGrp4 {c |}
{space 10}2  {c |}{col 14}{res}{space 2} .8300924{col 26}{space 2} .0227712{col 37}{space 1}   36.45{col 46}{space 3}0.000{col 54}{space 4} .7854016{col 67}{space 3} .8747832
{txt}{space 10}3  {c |}{col 14}{res}{space 2} .9510561{col 26}{space 2} .0280033{col 37}{space 1}   33.96{col 46}{space 3}0.000{col 54}{space 4} .8960955{col 67}{space 3} 1.006017
{txt}{space 10}4  {c |}{col 14}{res}{space 2}  .797854{col 26}{space 2} .0358343{col 37}{space 1}   22.27{col 46}{space 3}0.000{col 54}{space 4} .7275024{col 67}{space 3} .8682055
{txt}{space 12} {c |}
std_same_p~4 {c |}{col 14}{res}{space 2}-.0174547{col 26}{space 2}  .005989{col 37}{space 1}   -2.91{col 46}{space 3}0.004{col 54}{space 4}-.0292085{col 67}{space 3}-.0057009
{txt}{space 12} {c |}
{space 7}Race5 {c |}
{space 10}2  {c |}{col 14}{res}{space 2}-1.147644{col 26}{space 2} .0302037{col 37}{space 1}  -38.00{col 46}{space 3}0.000{col 54}{space 4}-1.206921{col 67}{space 3}-1.088366
{txt}{space 10}3  {c |}{col 14}{res}{space 2}-.1090002{col 26}{space 2} .0722856{col 37}{space 1}   -1.51{col 46}{space 3}0.132{col 54}{space 4}-.2508652{col 67}{space 3} .0328649
{txt}{space 10}4  {c |}{col 14}{res}{space 2}-.4053723{col 26}{space 2} .0519696{col 37}{space 1}   -7.80{col 46}{space 3}0.000{col 54}{space 4} -.507366{col 67}{space 3}-.3033786
{txt}{space 10}5  {c |}{col 14}{res}{space 2}-.7678096{col 26}{space 2} .0397165{col 37}{space 1}  -19.33{col 46}{space 3}0.000{col 54}{space 4}-.8457559{col 67}{space 3}-.6898634
{txt}{space 12} {c |}
std_same_~e5 {c |}{col 14}{res}{space 2}  .004527{col 26}{space 2} .0095081{col 37}{space 1}    0.48{col 46}{space 3}0.634{col 54}{space 4}-.0141333{col 67}{space 3} .0231872
{txt}{space 3}1.BornUSA {c |}{col 14}{res}{space 2} .4526259{col 26}{space 2} .0378773{col 37}{space 1}   11.95{col 46}{space 3}0.000{col 54}{space 4} .3782894{col 67}{space 3} .5269624
{txt}std_same_p~A {c |}{col 14}{res}{space 2} .0056886{col 26}{space 2} .0072761{col 37}{space 1}    0.78{col 46}{space 3}0.435{col 54}{space 4}-.0085913{col 67}{space 3} .0199685
{txt}{space 12} {c |}
{space 4}MarStat5 {c |}
{space 10}2  {c |}{col 14}{res}{space 2} .6836293{col 26}{space 2} .0215387{col 37}{space 1}   31.74{col 46}{space 3}0.000{col 54}{space 4} .6413567{col 67}{space 3} .7259019
{txt}{space 10}3  {c |}{col 14}{res}{space 2} .9596247{col 26}{space 2} .0143059{col 37}{space 1}   67.08{col 46}{space 3}0.000{col 54}{space 4} .9315353{col 67}{space 3}  .987714
{txt}{space 10}4  {c |}{col 14}{res}{space 2}-.3088691{col 26}{space 2} .0916434{col 37}{space 1}   -3.37{col 46}{space 3}0.001{col 54}{space 4}-.4887264{col 67}{space 3}-.1290118
{txt}{space 10}5  {c |}{col 14}{res}{space 2}  .636664{col 26}{space 2} .0187825{col 37}{space 1}   33.90{col 46}{space 3}0.000{col 54}{space 4} .5997957{col 67}{space 3} .6735323
{txt}{space 12} {c |}
std_same_~t5 {c |}{col 14}{res}{space 2}-.0255763{col 26}{space 2} .0058647{col 37}{space 1}   -4.36{col 46}{space 3}0.000{col 54}{space 4}-.0370863{col 67}{space 3}-.0140662
{txt}{space 1}Rat_Poverty {c |}{col 14}{res}{space 2}-.0047919{col 26}{space 2}  .001854{col 37}{space 1}   -2.58{col 46}{space 3}0.010{col 54}{space 4}-.0084305{col 67}{space 3}-.0011533
{txt}{space 1}Rat_Mig_Cum {c |}{col 14}{res}{space 2}-.7197979{col 26}{space 2}  .462591{col 37}{space 1}   -1.56{col 46}{space 3}0.120{col 54}{space 4}-1.627663{col 67}{space 3} .1880669
{txt}{space 5}Pop_Den {c |}{col 14}{res}{space 2}-1.88e-06{col 26}{space 2} 4.33e-06{col 37}{space 1}   -0.43{col 46}{space 3}0.664{col 54}{space 4}-.0000104{col 67}{space 3} 6.62e-06
{txt}{space 1}Rat_GC_ProE {c |}{col 14}{res}{space 2}-9.17e-07{col 26}{space 2} 8.79e-06{col 37}{space 1}   -0.10{col 46}{space 3}0.917{col 54}{space 4}-.0000182{col 67}{space 3} .0000163
{txt}{space 1}Rat_GC_ProM {c |}{col 14}{res}{space 2} .0000122{col 26}{space 2} .0000109{col 37}{space 1}    1.12{col 46}{space 3}0.264{col 54}{space 4}-9.24e-06{col 67}{space 3} .0000337
{txt}{space 1}Rat_GC_ProB {c |}{col 14}{res}{space 2}-2.77e-06{col 26}{space 2} .0000383{col 37}{space 1}   -0.07{col 46}{space 3}0.942{col 54}{space 4}-.0000779{col 67}{space 3} .0000723
{txt}{space 2}Rat_GC_Jew {c |}{col 14}{res}{space 2} .0001104{col 26}{space 2} .0000669{col 37}{space 1}    1.65{col 46}{space 3}0.099{col 54}{space 4}-.0000209{col 67}{space 3} .0002417
{txt}{space 2}Rat_GC_Oth {c |}{col 14}{res}{space 2}-.0000732{col 26}{space 2} .0000209{col 37}{space 1}   -3.50{col 46}{space 3}0.000{col 54}{space 4}-.0001142{col 67}{space 3}-.0000321
{txt}{space 4}1.UnEmpl {c |}{col 14}{res}{space 2}  1.92543{col 26}{space 2} .0443574{col 37}{space 1}   43.41{col 46}{space 3}0.000{col 54}{space 4} 1.837139{col 67}{space 3} 2.013721
{txt}std_same_p~l {c |}{col 14}{res}{space 2} .0070872{col 26}{space 2} .0165409{col 37}{space 1}    0.43{col 46}{space 3}0.669{col 54}{space 4}  -.02546{col 67}{space 3} .0396345
{txt}{space 2}1.PhysProb {c |}{col 14}{res}{space 2} .3416501{col 26}{space 2} .0247124{col 37}{space 1}   13.83{col 46}{space 3}0.000{col 54}{space 4} .2931368{col 67}{space 3} .3901635
{txt}std_same_p~b {c |}{col 14}{res}{space 2} -.127903{col 26}{space 2}  .012234{col 37}{space 1}  -10.45{col 46}{space 3}0.000{col 54}{space 4}-.1519207{col 67}{space 3}-.1038853
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-9.422722{col 26}{space 2} .1074588{col 37}{space 1}  -87.69{col 46}{space 3}0.000{col 54}{space 4} -9.63362{col 67}{space 3}-9.211823
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}{com}.         estimates store m1 
. 
.         if (`model_version' <= 4){c -(}
.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi Female "binary" "mi"

                           {txt}Female
{hline 61}
      Percentiles      Smallest
 1%    {res}        0              0
{txt} 5%    {res}        0              0
{txt}10%    {res}        0              0       {txt}Obs         {res}135,582,733
{txt}25%    {res}        0              0       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res}        1                      {txt}Mean          {res} .5158124
                        {txt}Largest       Std. Dev.     {res} .4997499
{txt}75%    {res}        1              1
{txt}90%    {res}        1              1       {txt}Variance      {res}   .24975
{txt}95%    {res}        1              1       {txt}Skewness      {res}-.0632814
{txt}99%    {res}        1              1       {txt}Kurtosis      {res} 1.004005

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0028
{txt}{col 47}Largest FMI{col 65}= {res}      0.0045
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      909.14
{txt}{col 47}        avg{col 65}= {res}      911.25
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      913.36

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:Female}{space 10}{txt:=} {space 10}0}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:Female}{space 10}{txt:=} {space 10}1}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0002386{col 26}{space 2} 1.89e-06{col 37}{space 1}  126.45{col 46}{space 3}0.000{col 54}{space 4} .0002348{col 67}{space 3} .0002423
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0000628{col 26}{space 2} 8.01e-07{col 37}{space 1}   78.39{col 46}{space 3}0.000{col 54}{space 4} .0000612{col 67}{space 3} .0000643
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0001758}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000021}}}{space 1}{space 1}{ralign 8:{res:{sf:-8.3e+01}}}{space 1}{space 1}{ralign 8:{res:{sf: 0.0e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0001799}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0001716}}}{space 1}
{com}.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi AgeGrp4 "categorical" "mi"
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Pairwise comparisons of predictive margins{col 47}Number of obs{col 65}= {res}  11,079,265

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0185
{txt}{col 47}Largest FMI{col 65}= {res}      0.0680
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      599.02
{txt}{col 47}        avg{col 65}= {res}      825.74
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      910.28

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 14}{hline 12}
{col 14}{c |}   Contrast{col 26}   Std. Err.{col 38}     [95% Con{col 51}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 14}{hline 12}
{space 5}AgeGrp4 {c |}
{space 5}2 vs 1  {c |}{col 14}{res}{space 2} .0000939{col 26}{space 2} 2.07e-06{col 37}{space 5} .0000898{col 51}{space 3} .0000979
{txt}{space 5}3 vs 1  {c |}{col 14}{res}{space 2} .0001152{col 26}{space 2} 2.94e-06{col 37}{space 5} .0001095{col 51}{space 3}  .000121
{txt}{space 5}4 vs 1  {c |}{col 14}{res}{space 2} .0000886{col 26}{space 2} 4.13e-06{col 37}{space 5} .0000805{col 51}{space 3} .0000967
{txt}{space 5}3 vs 2  {c |}{col 14}{res}{space 2} .0000214{col 26}{space 2} 2.47e-06{col 37}{space 5} .0000165{col 51}{space 3} .0000262
{txt}{space 5}4 vs 2  {c |}{col 14}{res}{space 2}-5.28e-06{col 26}{space 2} 3.92e-06{col 37}{space 5} -.000013{col 51}{space 3} 2.43e-06
{txt}{space 5}4 vs 3  {c |}{col 14}{res}{space 2}-.0000267{col 26}{space 2} 3.28e-06{col 37}{space 5}-.0000331{col 51}{space 3}-.0000202
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 14}{hline 12}
{res}{com}.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi Race5 "categorical" "mi"
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Pairwise comparisons of predictive margins{col 47}Number of obs{col 65}= {res}  11,079,265

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0024
{txt}{col 47}Largest FMI{col 65}= {res}      0.0050
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      908.09
{txt}{col 47}        avg{col 65}= {res}      911.88
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      914.12

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 14}{hline 12}
{col 14}{c |}   Contrast{col 26}   Std. Err.{col 38}     [95% Con{col 51}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 14}{hline 12}
{space 7}Race5 {c |}
{space 5}2 vs 1  {c |}{col 14}{res}{space 2}-.0001217{col 26}{space 2} 2.48e-06{col 37}{space 5}-.0001265{col 51}{space 3}-.0001168
{txt}{space 5}3 vs 1  {c |}{col 14}{res}{space 2}-.0000184{col 26}{space 2} .0000115{col 37}{space 5}-.0000411{col 51}{space 3} 4.25e-06
{txt}{space 5}4 vs 1  {c |}{col 14}{res}{space 2}-.0000594{col 26}{space 2} 6.44e-06{col 37}{space 5} -.000072{col 51}{space 3}-.0000468
{txt}{space 5}5 vs 1  {c |}{col 14}{res}{space 2}-.0000955{col 26}{space 2} 3.74e-06{col 37}{space 5}-.0001029{col 51}{space 3}-.0000882
{txt}{space 5}3 vs 2  {c |}{col 14}{res}{space 2} .0001033{col 26}{space 2}  .000012{col 37}{space 5} .0000797{col 51}{space 3} .0001268
{txt}{space 5}4 vs 2  {c |}{col 14}{res}{space 2} .0000623{col 26}{space 2} 5.49e-06{col 37}{space 5} .0000515{col 51}{space 3} .0000731
{txt}{space 5}5 vs 2  {c |}{col 14}{res}{space 2} .0000262{col 26}{space 2} 3.13e-06{col 37}{space 5}   .00002{col 51}{space 3} .0000323
{txt}{space 5}4 vs 3  {c |}{col 14}{res}{space 2} -.000041{col 26}{space 2}  .000014{col 37}{space 5}-.0000684{col 51}{space 3}-.0000136
{txt}{space 5}5 vs 3  {c |}{col 14}{res}{space 2}-.0000771{col 26}{space 2} .0000127{col 37}{space 5} -.000102{col 51}{space 3}-.0000522
{txt}{space 5}5 vs 4  {c |}{col 14}{res}{space 2}-.0000361{col 26}{space 2} 5.80e-06{col 37}{space 5}-.0000475{col 51}{space 3}-.0000248
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 14}{hline 12}
{res}{com}.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi BornUSA "binary" "mi"

                           {txt}BornUSA
{hline 61}
      Percentiles      Smallest
 1%    {res}        0              0
{txt} 5%    {res}        0              0
{txt}10%    {res}        1              0       {txt}Obs         {res}135,582,733
{txt}25%    {res}        1              0       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res}        1                      {txt}Mean          {res} .9265123
                        {txt}Largest       Std. Dev.     {res} .2609353
{txt}75%    {res}        1              1
{txt}90%    {res}        1              1       {txt}Variance      {res} .0680872
{txt}95%    {res}        1              1       {txt}Skewness      {res}-3.269105
{txt}99%    {res}        1              1       {txt}Kurtosis      {res} 11.68704

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0005
{txt}{col 47}Largest FMI{col 65}= {res}      0.0010
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      914.00
{txt}{col 47}        avg{col 65}= {res}      914.41
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      914.83

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:BornUSA}{space 9}{txt:=} {space 10}0}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:BornUSA}{space 9}{txt:=} {space 10}1}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0000974{col 26}{space 2} 3.37e-06{col 37}{space 1}   28.92{col 46}{space 3}0.000{col 54}{space 4} .0000908{col 67}{space 3}  .000104
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001532{col 26}{space 2} 1.13e-06{col 37}{space 1}  136.06{col 46}{space 3}0.000{col 54}{space 4}  .000151{col 67}{space 3} .0001554
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000557}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000039}}}{space 1}{space 1}{ralign 8:{res:{sf: 1.4e+01}}}{space 1}{space 1}{ralign 8:{res:{sf: 0.0e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000480}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000634}}}{space 1}
{com}.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi MarStat5 "categorical" "mi"
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Pairwise comparisons of predictive margins{col 47}Number of obs{col 65}= {res}  11,079,265

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0163
{txt}{col 47}Largest FMI{col 65}= {res}      0.0258
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      837.30
{txt}{col 47}        avg{col 65}= {res}      874.51
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      910.62

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 14}{hline 12}
{col 14}{c |}   Contrast{col 26}   Std. Err.{col 38}     [95% Con{col 51}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 14}{hline 12}
{space 4}MarStat5 {c |}
{space 5}2 vs 1  {c |}{col 14}{res}{space 2} .0000992{col 26}{space 2} 4.16e-06{col 37}{space 5} .0000911{col 51}{space 3} .0001074
{txt}{space 5}3 vs 1  {c |}{col 14}{res}{space 2} .0001629{col 26}{space 2} 3.10e-06{col 37}{space 5} .0001568{col 51}{space 3}  .000169
{txt}{space 5}4 vs 1  {c |}{col 14}{res}{space 2}-.0000269{col 26}{space 2} 6.90e-06{col 37}{space 5}-.0000404{col 51}{space 3}-.0000134
{txt}{space 5}5 vs 1  {c |}{col 14}{res}{space 2} .0000901{col 26}{space 2} 2.90e-06{col 37}{space 5} .0000844{col 51}{space 3} .0000958
{txt}{space 5}3 vs 2  {c |}{col 14}{res}{space 2} .0000637{col 26}{space 2} 4.91e-06{col 37}{space 5}  .000054{col 51}{space 3} .0000733
{txt}{space 5}4 vs 2  {c |}{col 14}{res}{space 2}-.0001261{col 26}{space 2} 8.02e-06{col 37}{space 5}-.0001419{col 51}{space 3}-.0001104
{txt}{space 5}5 vs 2  {c |}{col 14}{res}{space 2}-9.19e-06{col 26}{space 2} 5.22e-06{col 37}{space 5}-.0000194{col 51}{space 3} 1.06e-06
{txt}{space 5}4 vs 3  {c |}{col 14}{res}{space 2}-.0001898{col 26}{space 2} 7.01e-06{col 37}{space 5}-.0002036{col 51}{space 3} -.000176
{txt}{space 5}5 vs 3  {c |}{col 14}{res}{space 2}-.0000729{col 26}{space 2} 3.87e-06{col 37}{space 5}-.0000805{col 51}{space 3}-.0000653
{txt}{space 5}5 vs 4  {c |}{col 14}{res}{space 2} .0001169{col 26}{space 2} 6.66e-06{col 37}{space 5} .0001039{col 51}{space 3}   .00013
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 14}{hline 12}
{res}{com}.                 
.                 if (`model_version' == 1 | `model_version' == 2) {c -(}
.                         estimates restore m1
.                         mchange_mi RAT_Female "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_AgeGrp4_2 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_AgeGrp4_3 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_AgeGrp4_4 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_Race5_2 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_Race5_3 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_Race5_4 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_Race5_5 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_BornUSA "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_MarStat5_2 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_MarStat5_3 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_MarStat5_4 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_MarStat5_5 "continuous" "mi"
.                 {c )-}
.                 if (`model_version' == 3 | `model_version' == 4) {c -(}
.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi std_same_prop_Sex "continuous" "mi"

                      {txt}std_same_prop_Sex
{hline 61}
      Percentiles      Smallest
 1%    {res}-2.509072      -4.036727
{txt} 5%    {res}-1.338853      -4.036727
{txt}10%    {res}-.9852811      -4.036727       {txt}Obs         {res}135,582,733
{txt}25%    {res}-.5357422      -4.036727       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} .0312444                      {txt}Mean          {res}  .027466
                        {txt}Largest       Std. Dev.     {res} .9119966
{txt}75%    {res} .5763232       4.036727
{txt}90%    {res} 1.032018       4.036727       {txt}Variance      {res} .8317377
{txt}95%    {res} 1.401881       4.036727       {txt}Skewness      {res} .0293361
{txt}99%    {res} 2.682605       4.036727       {txt}Kurtosis      {res} 4.966479

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0006
{txt}{col 47}Largest FMI{col 65}= {res}      0.0012
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      913.79
{txt}{col 47}        avg{col 65}= {res}      913.83
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      913.88

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:std_same_p~x}{space 4}{txt:=} {space 2}-4.036727}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:std_same_p~x}{space 4}{txt:=} {space 3}4.036727}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001579{col 26}{space 2} 4.83e-06{col 37}{space 1}   32.69{col 46}{space 3}0.000{col 54}{space 4} .0001484{col 67}{space 3} .0001673
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001391{col 26}{space 2} 4.35e-06{col 37}{space 1}   32.01{col 46}{space 3}0.000{col 54}{space 4} .0001306{col 67}{space 3} .0001476
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0000188}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000090}}}{space 1}{space 1}{ralign 8:{res:{sf:-2.1e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0367338}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000364}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000012}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi std_same_prop_AgeGrp4 "continuous" "mi"

                    {txt}std_same_prop_AgeGrp4
{hline 61}
      Percentiles      Smallest
 1%    {res}-2.469153      -5.047376
{txt} 5%    {res}-1.445343      -5.047376
{txt}10%    {res}-1.068901      -5.047376       {txt}Obs         {res}135,582,733
{txt}25%    {res} -.504598      -5.047376       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} .0808892                      {txt}Mean          {res}  .127631
                        {txt}Largest       Std. Dev.     {res} 1.053953
{txt}75%    {res} .6954771       6.771116
{txt}90%    {res} 1.347035       6.771116       {txt}Variance      {res} 1.110817
{txt}95%    {res} 1.891168       6.771116       {txt}Skewness      {res} .3782019
{txt}99%    {res} 2.953919       6.771116       {txt}Kurtosis      {res}  5.03041

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0008
{txt}{col 47}Largest FMI{col 65}= {res}      0.0016
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      913.27
{txt}{col 47}        avg{col 65}= {res}      913.36
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      913.45

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:std_same_p~4}{space 4}{txt:=} {space 2}-5.047376}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:std_same_p~4}{space 4}{txt:=} {space 3}6.771116}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001626{col 26}{space 2} 5.02e-06{col 37}{space 1}   32.39{col 46}{space 3}0.000{col 54}{space 4} .0001527{col 67}{space 3} .0001724
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001323{col 26}{space 2} 5.43e-06{col 37}{space 1}   24.37{col 46}{space 3}0.000{col 54}{space 4} .0001216{col 67}{space 3} .0001429
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0000303}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000103}}}{space 1}{space 1}{ralign 8:{res:{sf:-2.9e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0032815}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000504}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000101}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi std_same_prop_Race5 "continuous" "mi"

                     {txt}std_same_prop_Race5
{hline 61}
      Percentiles      Smallest
 1%    {res}-2.270736      -3.614547
{txt} 5%    {res}-1.387242      -3.614547
{txt}10%    {res}-.9511417      -3.614547       {txt}Obs         {res}135,582,733
{txt}25%    {res}-.2595515      -3.614547       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} .4283832                      {txt}Mean          {res} .4242603
                        {txt}Largest       Std. Dev.     {res} 1.233905
{txt}75%    {res} 1.013094       8.750521
{txt}90%    {res} 1.238259       8.750521       {txt}Variance      {res} 1.522523
{txt}95%    {res} 2.074838       8.750521       {txt}Skewness      {res} 1.661957
{txt}99%    {res} 6.478749       8.750521       {txt}Kurtosis      {res} 11.30467

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0013
{txt}{col 47}Largest FMI{col 65}= {res}      0.0026
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      911.97
{txt}{col 47}        avg{col 65}= {res}      911.97
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      911.98

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:std_same_~e5}{space 4}{txt:=} {space 2}-3.614547}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:std_same_~e5}{space 4}{txt:=} {space 3}8.750521}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001459{col 26}{space 2} 5.28e-06{col 37}{space 1}   27.65{col 46}{space 3}0.000{col 54}{space 4} .0001356{col 67}{space 3} .0001563
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001543{col 26}{space 2} .0000127{col 37}{space 1}   12.17{col 46}{space 3}0.000{col 54}{space 4} .0001294{col 67}{space 3} .0001792
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000084}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000178}}}{space 1}{space 1}{ralign 8:{res:{sf:0.4708513}}}{space 1}{space 1}{ralign 8:{res:{sf:0.6378589}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000266}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000434}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi std_same_prop_BornUSA "continuous" "mi"

                    {txt}std_same_prop_BornUSA
{hline 61}
      Percentiles      Smallest
 1%    {res}-5.045306      -7.936538
{txt} 5%    {res}-2.723737      -7.936538
{txt}10%    {res}-1.754185      -7.936538       {txt}Obs         {res}135,582,733
{txt}25%    {res}-.5298624      -7.936538       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} .1534955                      {txt}Mean          {res}-.1095675
                        {txt}Largest       Std. Dev.     {res} 1.400314
{txt}75%    {res} .6318931       7.936538
{txt}90%    {res} .8333296       7.936538       {txt}Variance      {res} 1.960879
{txt}95%    {res} .9703126       7.936538       {txt}Skewness      {res}-.7379761
{txt}99%    {res} 4.547946       7.936538       {txt}Kurtosis      {res} 8.840372

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0007
{txt}{col 47}Largest FMI{col 65}= {res}      0.0014
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      913.53
{txt}{col 47}        avg{col 65}= {res}      913.54
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      913.55

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:std_same_p~A}{space 4}{txt:=} {space 2}-7.936538}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:std_same_p~A}{space 4}{txt:=} {space 3}7.936538}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001424{col 26}{space 2} 7.72e-06{col 37}{space 1}   18.45{col 46}{space 3}0.000{col 54}{space 4} .0001272{col 67}{space 3} .0001575
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001558{col 26}{space 2} 9.65e-06{col 37}{space 1}   16.15{col 46}{space 3}0.000{col 54}{space 4} .0001369{col 67}{space 3} .0001747
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000134}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000173}}}{space 1}{space 1}{ralign 8:{res:{sf:0.7789207}}}{space 1}{space 1}{ralign 8:{res:{sf:0.4362274}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000204}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000473}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi std_same_prop_MarStat5 "continuous" "mi"

                   {txt}std_same_prop_MarStat5
{hline 61}
      Percentiles      Smallest
 1%    {res}-2.292691      -4.824302
{txt} 5%    {res}-1.574807      -4.824302
{txt}10%    {res}-1.168543      -4.824302       {txt}Obs         {res}135,582,272
{txt}25%    {res}-.5012236      -4.824302       {txt}Sum of Wgt. {res}  135582272

{txt}50%    {res} .1833285                      {txt}Mean          {res} .1423099
                        {txt}Largest       Std. Dev.     {res} 1.051846
{txt}75%    {res} .8066256         5.1647
{txt}90%    {res} 1.365528         5.1647       {txt}Variance      {res}  1.10638
{txt}95%    {res} 1.765109         5.1647       {txt}Skewness      {res} .0184338
{txt}99%    {res} 2.683285         5.1647       {txt}Kurtosis      {res}  4.57374

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0038
{txt}{col 47}Largest FMI{col 65}= {res}      0.0073
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      903.45
{txt}{col 47}        avg{col 65}= {res}      903.70
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      903.95

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:std_same_~t5}{space 4}{txt:=} {space 2}-4.824302}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:std_same_~t5}{space 4}{txt:=} {space 5}5.1647}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001686{col 26}{space 2} 5.08e-06{col 37}{space 1}   33.18{col 46}{space 3}0.000{col 54}{space 4} .0001586{col 67}{space 3} .0001786
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001306{col 26}{space 2} 3.89e-06{col 37}{space 1}   33.54{col 46}{space 3}0.000{col 54}{space 4}  .000123{col 67}{space 3} .0001383
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0000380}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000088}}}{space 1}{space 1}{ralign 8:{res:{sf:-4.3e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000164}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000552}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000208}}}{space 1}
{com}.                 {c )-}
.         
.                 if (`model_version' == 2 | `model_version' == 4){c -(}
.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi UnEmpl "categorical" "mi"
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Pairwise comparisons of predictive margins{col 47}Number of obs{col 65}= {res}  11,079,265

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.2377
{txt}{col 47}Largest FMI{col 65}= {res}      0.3336
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}       77.80
{txt}{col 47}        avg{col 65}= {res}      118.98
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      160.17

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 14}{hline 12}
{col 14}{c |}   Contrast{col 26}   Std. Err.{col 38}     [95% Con{col 51}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 14}{hline 12}
{space 6}UnEmpl {c |}
{space 5}1 vs 0  {c |}{col 14}{res}{space 2} .0006714{col 26}{space 2} .0000275{col 37}{space 5} .0006165{col 51}{space 3} .0007262
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 14}{hline 12}
{res}{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi PhysProb "categorical" "mi"
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Pairwise comparisons of predictive margins{col 47}Number of obs{col 65}= {res}  11,079,265

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0216
{txt}{col 47}Largest FMI{col 65}= {res}      0.0379
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      774.21
{txt}{col 47}        avg{col 65}= {res}      827.20
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      880.19

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 14}{hline 12}
{col 14}{c |}   Contrast{col 26}   Std. Err.{col 38}     [95% Con{col 51}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 14}{hline 12}
{space 4}PhysProb {c |}
{space 5}1 vs 0  {c |}{col 14}{res}{space 2} .0000565{col 26}{space 2} 4.55e-06{col 37}{space 5} .0000476{col 51}{space 3} .0000655
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 14}{hline 12}
{res}{com}.                 {c )-}
.         
.                 if (`model_version' == 2){c -(}
.                         estimates restore m1
.                         mchange_mi RAT_UnEmpl "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_PhysProb "continuous" "mi"
.                 {c )-}
.                 if (`model_version' == 4){c -(}
.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi std_same_prop_UnEmpl "continuous" "mi"

                    {txt}std_same_prop_UnEmpl
{hline 61}
      Percentiles      Smallest
 1%    {res}-2.422866      -5.256251
{txt} 5%    {res}-1.473278      -5.256251
{txt}10%    {res}-1.055818      -5.256251       {txt}Obs         {res}135,334,296
{txt}25%    {res}-.6145594      -5.256251       {txt}Sum of Wgt. {res}  135334296

{txt}50%    {res} .1522239                      {txt}Mean          {res} .0244813
                        {txt}Largest       Std. Dev.     {res} .9458323
{txt}75%    {res} .6237007       5.256251
{txt}90%    {res} 1.075561       5.256251       {txt}Variance      {res} .8945988
{txt}95%    {res} 1.334634       5.256251       {txt}Skewness      {res}-1.054592
{txt}99%    {res} 1.864986       5.256251       {txt}Kurtosis      {res}  7.51882

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0825
{txt}{col 47}Largest FMI{col 65}= {res}      0.1378
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      306.73
{txt}{col 47}        avg{col 65}= {res}      310.60
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      314.46

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:std_same_p~l}{space 4}{txt:=} {space 2}-5.256251}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:std_same_p~l}{space 4}{txt:=} {space 3}5.256251}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001431{col 26}{space 2} .0000125{col 37}{space 1}   11.45{col 46}{space 3}0.000{col 54}{space 4} .0001185{col 67}{space 3} .0001677
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001541{col 26}{space 2} .0000134{col 37}{space 1}   11.50{col 46}{space 3}0.000{col 54}{space 4} .0001278{col 67}{space 3} .0001805
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000111}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000258}}}{space 1}{space 1}{ralign 8:{res:{sf:0.4282625}}}{space 1}{space 1}{ralign 8:{res:{sf:0.6685606}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000396}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000617}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi std_same_prop_PhysProb "continuous" "mi"

                   {txt}std_same_prop_PhysProb
{hline 61}
      Percentiles      Smallest
 1%    {res} -2.23507      -4.051959
{txt} 5%    {res}-1.302234      -4.051959
{txt}10%    {res}-.9170214      -4.051959       {txt}Obs         {res}135,576,336
{txt}25%    {res}-.4860743      -4.051959       {txt}Sum of Wgt. {res}  135576336

{txt}50%    {res} .3726953                      {txt}Mean          {res} .2796852
                        {txt}Largest       Std. Dev.     {res} .9714683
{txt}75%    {res} .9583569       4.051959
{txt}90%    {res}  1.46031       4.051959       {txt}Variance      {res} .9437506
{txt}95%    {res} 1.737546       4.051959       {txt}Skewness      {res}-.3642149
{txt}99%    {res} 2.235158       4.051959       {txt}Kurtosis      {res} 3.305471

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0230
{txt}{col 47}Largest FMI{col 65}= {res}      0.0442
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      737.79
{txt}{col 47}        avg{col 65}= {res}      742.63
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      747.47

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:std_same_p~b}{space 4}{txt:=} {space 2}-4.051959}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:std_same_p~b}{space 4}{txt:=} {space 3}4.051959}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0002627{col 26}{space 2} .0000143{col 37}{space 1}   18.42{col 46}{space 3}0.000{col 54}{space 4} .0002347{col 67}{space 3} .0002906
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0000932{col 26}{space 2} 4.25e-06{col 37}{space 1}   21.91{col 46}{space 3}0.000{col 54}{space 4} .0000849{col 67}{space 3} .0001016
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0001694}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000184}}}{space 1}{space 1}{ralign 8:{res:{sf:-9.2e+00}}}{space 1}{space 1}{ralign 8:{res:{sf: 0.0e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0002055}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0001333}}}{space 1}
{com}.                 {c )-}               
.         {c )-}
. {c )-}
{txt}
{com}. 
. if ("`model'" == "logit"){c -(}
.         svy: logit Suic i.St `model_eq', or 
.         estimates store m1 
. 
.         if (`model_version' <= 4){c -(}
.                 mchange `margin_demographics', amount(all) delta(100) statistics(all) decimals(7)
.         {c )-}
. 
. {c )-}
{txt}
{com}. 
. if (`model_version' == 5 | `model_version' == 6) {c -(}
.         estimates restore m1
.         margin_interact std_same_prop_Sex Female 5 `model'
.         estimates restore m1
.         margin_interact std_same_prop_AgeGrp4 AgeGrp4 5 `model'
.         estimates restore m1
.         margin_interact std_same_prop_Race5 Race5 5 `model'
.         estimates restore m1
.         margin_interact std_same_prop_BornUSA BornUSA 5 `model'
.         estimates restore m1
.         margin_interact std_same_prop_MarStat5 MarStat5 5 `model'
. 
.         if (`model_version' == 6) {c -(}
.                 estimates restore m1
.                 margin_interact std_same_prop_UnEmpl UnEmpl 5 `model'
.                 estimates restore m1
.                 margin_interact std_same_prop_PhysProb PhysProb 5 `model'
.         {c )-}
. 
. {c )-}
{txt}
{com}. 
. 
. log close 
      {txt}name:  {res}<unnamed>
       {txt}log:  {res}/N/project/suicide_study/pnas_replication/results/log/mi_4.smcl
  {txt}log type:  {res}smcl
 {txt}closed on:  {res}20 Aug 2020, 23:00:26
{txt}{.-}
{smcl}
{txt}{sf}{ul off}